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ABSTRACT 


Lunisolar  perturbations  for  general  terms  of  the  disturbing  function  were 
derived  by  Kaula  (1962).  However,  his  formulas  use  equatorial  elements  for  the  Moon 
and  do  not  give  a  definite  algorithm  for  computational  procedures.  As  Kozai  (1966) 
suggested,  both  inclination  and  node  of  the  Moon's  orbit  with  respect  to  the  equator  of 
the  Earth  are  not  simple  functions  of  time,  while  the  same  elements  with  respect  to 
the  ecliptic  are  well  approximated  by  a  constant  and  a  linear  function  of  time,  res¬ 
pectively.  In  the  present  work,  we  obtain  the  disturbing  function  for  the  Moon’s  per¬ 
turbations  using  ecliptic  elements  for  the  Moon  and  equatorial  elements  for  the  satellite. 
Secular,  long-period,  and  short -period  perturbations  are  then  computed,  with  the 
expressions  kept  in  closed  form  in  both  inclination  and  eccentricity  of  the  satellite. 
Alternative  expressions  for  short-period  perturbations  of  high  satellites  are  also 
given,  assuming  small  values  of  the  eccentricity.  The  Moon's  position  is  specified 
by  the  inclination,  node,  argument  of  perigee,  true  (or  mean)  longitude,  and  its  radius 
vector  from  the  center  of  the  Earth.  We  can  then  apply  the  results  to  numerical 
integration  by  using  coordinates  of  the  Moon  from  ephemeris  tapes  or  to  analytical 
representation  by  using  results  from  lunar  theory,  with  the  Moon's  motion  represented 
by  a  processing  and  rotating  elliptical  orbit. 
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RESUME 


Kaula  (1962)  a  dgduit  les  perturbations  lunisolaires  pour  les  termes  ggngraux 
de  la  fonction  de  perturbation.  Toutefois,  ses  formules  emploient  des  dldments 
dquatoriaux  pour  la  lune  et  ne  donnent  pas  un  algorithme  dgfini  pour  les 
proc£d£s  de  calcul.  Comme  l'a  sugg^rd  Kozai  (1966),  1 'inclinaison  et  le  noeud 
de  l'orbite  lunaire  ne  sont  pas  de  simples  fonctions  du  temps  par  rapport  3 
l'orbite  terrestre,  tandis  que  ces  monies  dldments  peuvent  @tre  assimilds, 
respectivement,  3  une  constante  et  cl  une  fonction  lindaire  du  temps  par  rapport 
3  1 '^cliptique.  Dans  cette  ytude  nous  avons  obtenu  la  fonction  de  perturbation 
pour  les  perturbations  lunaires  en  employant  des  dldmerts  dcliptiques  pour 
la  lune  et  des  dldments  dquatoriaux  pour  le  satellite.  Nous  avons  alors  cal¬ 
culi  les  perturbations  sdculaires  et  les  perturbations  £  longues  et  courtes 
pdriodes,  les  expressions  de  1 'inclinaison  et  de  1 'eccentricity  du  satellite 
restant  ddlimitdes.  Nous  donnons  aussi  des  expressions  alternatives  pour  des 
perturbations  a  courtes  pdriodes  de  hauts  .atellites,  en  supposant  que  les 
valeurs  de  1 'eccentricity  sont  petites.  La  position  de  la  lune  est  donn^e  par 
1 ' inclinaison,  le  noeud,  l'argument  du  pyrigye,  la  longitude  vraie  (  ou  moyenne), 
et  son  vecteur  cl  partir  du  centre  de  la  terre.  Nous  pouvons  ensuite  appliquer 
les  r£sultats  soit  h  une  integration  numyrique,  en  nous  servant  des  coordonnyes 
de  la  lune  donnyes  pa*'-  les  yphymyrides,  soit  a  une  reprysentation  analytique 
en  nous  servant  des  rysultats  de  la  thyorie  lunaire,  le  mcuvement  de  la  lune 
ytant  reprysenty  par  une  orbite  elliptique  rotatoire  avec  prycession. 


KOHCnEKl 


SHayeHiiH  jiyHHo-cojiHeHHUx  B03Myi»eKit:t  ehh  oGiuiix  hjic-hob  bogm;/!.  a- 
wme:l  ffyHKiimi  6mjtm  sHseneHM  Kayna  (1962).  OnnaKo  b  ero  TopMy.nax 
i!cnoab3yiOTcH  3KBaTopii3JibHM6  ojieweHTbi  hjif  jivhm  it  ne  hootch  onre-z.0- 
JieHHOro  an ropuTua  hjir  npOBeneHitn  BbimtcjioHHii .  CornacMO  yKa3air/fi:.: 
Kosait  (1966),  Kan  HaKjicn e;i it o  ,  tek  ii  yaen  Jiynnoii  opThth  cthoci't oji  bh c 
seMHoro  aKBaropa  ne  kbjifiiotcf  npocTMMK  Ip/hkhhhmh  BpeweHn,  b  to  btomh 
KaK  cTHOCHTenbHO  annum  hki:  to  :*e  caubie  ojicmphtli  xopoiso  annpOKcn.Ht- 
pOBaHb:  COOTBeTCTBOiiHO  KOHCTaHTHOli  W  JlMHeKHoT:  (}\y H K U. V.  & Ii  BpeMOHIT. 

B  naHHoil  pa6oxe  mm  nonyuaeM  BOPHymawinyio  tfiyuKiimo  juih  jiyHHKX  eosy;/- 
meHX”,  I!CnOJ!b3yH  OnSMeHTbl  3  KJIMnTItKM  HJI>i  JiyHH  I!  3  KBaTOpiTaJlbHHe 
ojiaMeHTw  hjih  cnyTHHKa.  3aT6M  EbimrcnHKiTCH  bokobuo,  noni obpc-mohhho 
i:  KpaTKOBpeweHHHo  F.03i.:yae»-:>!H,  npiweM  ii  s  HaKncHenitii  i*  b  bkcuoftt  : - 
ciTTeTe  cnyTHiina  BMpaateHwa  nant:  r  biiup  saM any tmx  1or:.?yn.  B  paaore 
Tanxe  npMBeneHbi  anbTer hhthbhko  BhipaacemtH  iijih  k p a  t  k  o  b  p  o u a h h w x 
BOOMyiBOHIil?  enyTHKKOB,  FaXOflFXHXC;!  Ha  BHCOKO.i  OpOHTO,  npHHeM 
cneJiaHH  jionyxeHHH  o  tom ,  hto  3kcuohtphcmtct  v m cot  wajiue  3iraMeH  :>• . 
flojio^eHHe  JiyHH  onpenenneTcn  HaKJiOHemreM ,  ,v3jiom ,  arrywenTOM  nepitreo  , 
iicthhhoI?  (wjiii  crorHoiO  nonroTC.  it  pannycoM-BOKTor or  ji.vhm  ot  u^htto 
3ewnit.  3a.Tew  nojiywenHue  pesynmaTM  mh  mot  pm  npiiMenwTb  a  jib  nrcr.en- 
Horo  iiht e rpxpo b&hhh ,  ncnojibsya  KoopramaTbi  .n y h h  it 3  3anncoi'i  ojiewepuji 
na  neHTax,  t*n it  nun  aHamiTunecKorc  MonemtpcBaHitn ,  ncnojibayn  I'oaynb- 

TBTbl  Jl.VMHOH  TOOpi'lT,  nplIM  CM  JiBHKCHKr  JS'/HK  HpOUCTaBJl  CHO  B  B11JIO 

npeucccitpyioutoii  it  Brarammolicf?  snnnnTH’iPCKoli  opOr.TH . 


LUNAR  PERTURBATIONS  ON 
ARTIFICIAL  SATELLITES  OF  THE  EARTH 


Giorgio  E.  O.  Giacaglia 

1.  ELEMENTS  FOR  THE  MOON  AND  OTHER  QUANTITIES 

Let  T  be  the  time  in  centuries  of  36525  ephemeris  days  from  J.  D.  2415020. 0. 
The  following  values  will  be  adopted: 

Eccentricity  of  the  Moon: 

e  =  0. 054900489  . 

Inclination  of  the  M<  on: 

I(  =  5°  8 '43V  427 

sin  (1  /2)  =  0.  044886967  . 

Mass  ratio,  Moon  to  Earth: 

rnf/in„  =  0.0123001  • 

Mean  equatorial  parallax  of  the  Moon: 

=  57'2V70  , 

where 

=  arc  sin  (a@/a  )  . 

Mean  equatorial  radius  of  the  Earth: 

a  =  6378160  m 
e 

a^  =  perturbed  semimajor  axis  of  the  Moon's  orbit. 

Thi3  work  was  supported  in  part  by  contract  N00014-7  l-A-0110-0004  from  the  Office 
of  Naval  Research,  Mathematics  Program,  and  in  part  by  grant  NGR  09-015-002 
from  the  National  Aeronautics  and  Space  Administration. 
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Mean  motion  in  longitude  of  the  Moon: 

n^  =  13°10'34V889902  day-1  . 

Mean  anomaly  of  the  Moon*. 

M^  =  36°55'16V80  +  1724878768V03  T  +  25V61  1 2  +  0V0438  T3  . 

Argument  of  Moon's  perigee  from  ecliptic  node: 

=  -25°40’13V6G  +  14648522'.'51  T  -  37V 17  T2  -  0V0450  T3  . 

Ecliptic  longitude  of  Moon's  ascending  node: 

O  =  259°10'59V79  -  6962911V23  T  +  7V48  T2  +  0V0080  T3  . 

Obliquity  of  the  ecliptic: 

£  =  23°27'08V26  -  46V 845  T  -  0V0059  T2  +  0V00181  T3  . 

Explicit  and  precise  variations  of  £  in  terms  of  the  elements  of  the  Moon  and  the  Sun 
are  also  available: 

£  —  Eq  +  £2  +  dcj 

(e.g. ,  Connaissance  des  Temi  s,  1971). 


2.  THE  DISTURBING  FUNCTION 


The  disturbing  function  can  be  written  as 


Gm 
R  =  - i 


{COS  ifr')  , 

where  if>'  is  the  geocentric  elongation  of  the  satellite  from  the  Moon. 


(1) 


The  following  size  considerations  apply: 


m. 


Gm.  = - — 

l  mr  + 


2  3 


2  3 


n,  a,  =*  0. 0123  n„  a 


£  £ 


£  c 


also  written  as 


Gm  =  N2  a?  ,  N2  =*  1.  59  X  10~5  rev2  day"2  . 

£  £  £  m 

The  satellite  Keplerian  negative  energy  is 
F0  =  n2  a2/2  , 

so  that  the  relative  size  of  the  perturbing  force  function  is  given  by 
v=R/F0=  2N2/n2  . 

For  low  satellites  (T  =  90  min),  v  ^  1.  2  X  10  For  high  satellites  (T  =  24  h), 
v  =*  3. 18  X  10  .  It  follows  that,  in  the  above  range  of  periods,  for  moderate  eccen¬ 

tricities,  the  dominant  part  of  the  disturbing  function  of  a  satellite  is  due  to  the  Earth 
oblateness  (J2),  and  lunar  perturbations  are  about  second  order  with  respect  to  this. 


Let  a,  6  and  a',  6'  be  the  right  ascension  and  declination  of  the  satellite  and  of 
the  Moon,  respectively  (in  an  equatorial  system).  Jt  follows  that 


cos  if>r  =  cos  6  cos  6'  cos  (a -a')  +  sin  6  sin  6' 


Therefore, 


R=I>t  \‘Rt  ■ 


where 


or,  making  use  of  Legendre's  addition  theorem, 


.  |  +  1 

R  =  a^  (- M-M  V"*  €  "f.  ■"  P™  (sin  6')  p"1  (sin  6)  cos  m(a  -a')  , 

i  \a>  \r  /  Z_J  m  (i  +  m)!  f  '  ’  t  ' 


where 


( l  ,  m  =  0 
^2  ,  m  it  0 


A™  =  (-A  €  P™  (sin  6') 

f  \rJ  (jf  +  m)!  m  i  ' 


cos  ma  , 


B™  =  (-^)  ^  f  P™  (sin  &')  sin  ma'  . 

I  \r  /  (i  +  m)!  mi'  ' 
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We  can  write 


m=0 


Let 


v  =  u  +  i  =  argument  of  latitude  of  satellite, 

£2  =  longitude  (equatorial)  of  node  of  satellite, 
1=  inclination  of  satellite  to  Earth  equator. 

By  considering  the  relations 


cos  (a  -  f2)  cos  6  =  cos  v  , 
sin  (a  -  fi)  cos  6  =  sin  v  cos  I  , 
sin  6  =  sin  v  sin  I  , 

it  follows  that 


(6) 


(■) 


V“'(l f  £  £  FimP 

m=0  p=0 


o 


,  m 


m 

L-Bj?  - 


i  -  m  even 


i  -  m  odd 


co^  [(i  -  2p)  v  +  m£2]  + 


rRrm 

Am 
L  l  . 


i  -  m  even 


sin  ((j?  -  2p)  v  +  mS2]' 


(8) 


i  -  m  odd 


where  (Kaula,  1961) 


F  (D-V 

i  mp' '  i!  (i  -i)I  (£ 


„2i-2f  . 

2  .  f  -m-2i  T 

- tttt,  sin  I  X 

-  m-  2i)I 


x£(r)-sii£(i'mi;2l+i)(p-i'-ik)<-i>k'q 


(9) 


and 


q  =  [■  ~m],  the  integral  part  of  (I  -  m)/2, 

i  =  0,1,2,...,  min  (p,q), 
j  =  0, 1, 2, . . . ,  m, 

k  =  all  values  for  which  the  coefficient  is  not  zero;  that  is,  p  -  i  >  k. 


3.  ROTATION  OF  SPHERICAL  HARMONICS  FOR  THE  MOON 


The  spherical  harmonics  P^11  (sin  6')  eima  are  expressed  in  terms  of  6^  and  a^, 


the  Moon's  ecliptic  latitude  and  longitude. 


The  relations  are  given  by 


cos  6'  ela  =  cos  5^  cos  +  i  (cos  6^,  sin  cos  £  -  sin  b^  sin  £)  , 


C - 1 


sin  b'  =  cos  6^  sin  sin  £  +  sin  6^  cos  £ 


From  the  well-known  properties  of  spherical  harmonics  under  rotation,  we  can  write 


,  « — *  ira. 

„m  ,  .  l ma7  \  1  m, r  „r  ,  .  ,  .  C  , 

Pf  (sin  6')  e  =  2^  ai  <Sin  V 


in  r 

where  a  »  is  a  function  of  £  only,  <nd 


>;m« -<-*>“$*$  if* 


Using  the  orthogonality  conditions  of  spherical  harmonics,  we  have 


m, r  frgj^>  (i  -r)!  +f 
a(  4  it  (t  +  r)  l  J 


cos  b.  d6.  X 
t  ( 


X  J  P™(sin  6')  eimG  p£  (sin  6^)  e  ^  da,  . 


Introducing  Equations  (10)  and  computing  the  above  integral,  we  find  (see  also  Lee, 
1971) 


7 


9 


mtr._  (l  -r)I  i(m-r)ir/2  m,r 
at  (l-m)!  e  ui 


(11) 


where,  for  m  +  r  >0, 


TTm,  r  .  -  m  f  t  +  m\  (  £\ 

V  =<-1)  (/-r)\°os2/ 

'H  +r. 


m  +  r 


£\r_m 

sin  f)  X 


X  Fl-f  +r,  l  +r  +  1,  m  +  r  + 1;  cos 


and,  for  m  +  r  ^  0, 


TTm,r  ..J-r/i-mW  E\~m'r  /  .  £\m_r 

V  '+D  ^r)vc0s2j  r>5j 


2  F 

X  Fl-f  -r,  i  -  r+  1,  -m-r+  1;  cos 


X 


(12a) 


(12b) 


In  the  above  relations,  F  is  the  usual  hypergeometric  series  ,F  ,  defined  by 

ST  (a)n  (b)n  xn 
F(a,b,c;x)=^-^fr  , 


n=0 


where 


(a)n -  a(a+  1)  (a+  2)  ...  (a  +  n -  1)  , 

(a)0  =  1  • 

In  both  cases,  U™'  r  is  a  polynomial  in  sin  £/2,  cos  £/2  since  at  least  one  of  the 

parameters  a,b  is  negative,  and  the  above  series  terminate.  The  distinction  of  the 

cares  m  +  r  £  C  is  necessary  to  avoid  a  singularity  in  F  due  to  a  negative  value  of  ?. 

in  i* 

Considering  this  fact,  another  possible  form  for  ’  ,  valid  in  any  case,  is  found 
to  be 


8 


.  ..iii-i  /  ~>m  +  r  / 

(stoD 


d*  r  i  “  rn  . "t*  m .  / 1 o\ 

-jTT  [z  (z_1)  1  »  (13) 

dz 


where  z  =  cos  (£/2). 


Now,  let 


2  A™’r  =  U™,r  +  (-l)r  U ]"»  "r 


2  B™» r  -  U™» r  -  (-l)r  U™*  "r 


It  follows  that,  for  tn  even, 


1  +  1  m  f 


(f  +  m)! 


f  +  1  /_  nm  * 

Jim  =  &)  1(W  L  (!  -r):  *r  Bim,r  PI  <sin  V  sln  [r(“c  +  D]  ’ 

i  r=l 

(15a) 


and,  for  m  odd, 


/aA£+1  <-l)m+1< 

Arn  J-l]  — - 

t  \TrJ  (f  +  m)! 


*  +  1  .  .  m  f 

®T  -  &)  uTsr  Z ,£ - r)!  e  r A?’ r  pTi  <sin  V 008  [r(“< 
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Making  use  of  relations  (7)  for  the  Moon,  we  find  that,  for  m  even: 


'"-0 


+  1  (-1)"^  J- 


W  L  V  • 

r=  0  P=0 


c  +  l  ,  ,.m  £ 

la  Y 


i  YrJ  (£  +  m) 


"E''-r)",>rrZF<rpV8inV  ’ 

r=0  p=0 


and,  for  m  odd: 


£  + 1  m+ 1  £  £ 

“73 - r~ 1  ^  B™,rV'  F.  (T .)  sin  0J  , 

jf  VrJ  (f  +m)!  L  j  '  '  r  £  /  ^  f-'p  '  ^pr  ’ 


i  +  1  ,  ,,m  i 

ar\  (-1)  c  __ 
(T  \  '  m 


risyr  E  - r):  ‘r  Aim’ r  E  F»rP  <V  c“  V  > 

r=0  p=0 


where 


V=('‘2p)vc  +  r(nit  2' 


and  the  functions  (Ir)  are  defined  by  Equation  (9). 


Finally,  the  function  R  can  be  written,  with  Equations  \8),  (14),  (16),  and  (17) 


taken  into  account: 


.  1  /a  \*  +  1  *  *  *  *  (-1)™  e  t  (£  -  s)l 

=  J  fit  (-0  V  V  V  V  _ _ SL  s _ _ 


R<  =  a  it)  Ir 


EEEE 

m=0  s=0  p=0  q=0 


(£  +  m) 


Ffmp  (I)  Fisq 


y/  f",  i  rn  s  TTm«  -s  ,a  ,  a ^  s  /a  a ^  3 

X  [<■*>  V  cos<Vn  +  V)  +  V  008  <%m  -  V>J  > 


fpm  l qs  '  I pm  ^qs'. 


where 


9 .  =  (£  -  2p)  v  +  mi2  . 

ipm  ' 
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4.  SECULAR  AND  LONG-PERIOD  TERMS  OF  THE  DISTURBING  FUNCTION 


If  we  assume  that  no  resonance  occurs  between  the  orbital  motion  of  the  satellite 
and  that  of  the  Moon  -  that  is,  pn  +  p'n,  is  not  small  for  small  intergers  p,p'  not 
simultaneously  zero  -  then  the  elimination  of  short-period  terms  (depending  on  the 
mean  anomaly  of  the  satellite,  M)  from  the  disturbing  function  can.  be  obtained  by 
making  use  of  the  well-known  integrals 


2tt 

2 h  I  (at  sin(i-2p)fdM=0  , 


0 

2ir 


2 \\  ©'cos^rtfdMMl  +  pV*'1  <;'0‘2P(»=(^pV  "  Hfp(2p.{)(P) 


(21) 


In  the  above  relations, 


=  e^l  +  ^1  -e2^ 


(22) 


and  the  X’s  are  Hansen's  coefficients  (e.g.,  see  Plummer,  1960,  p.  45)  and  the  H's  are 
Kaula's  coefficients  (Kaula,  1961).  They  are  defined  by,  for  2p-i  >  0, 

Hip(2p-I) 

; 

: 

1 

and,  for  0, 


=  (-P) 


2p  -l 


/2p  + 
\2p  - 


A 

V 


FM  -1,  2p-2i  -1,  2p-*;  p  ) 


(23) 
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-i. ...  . 


In  both  cases,  they  are  polynomials  in  p.  The  distinction  again  is  necessary  in  order 
to  avoid  singular  representation. 


The  long-period  and  secular  part  of  the  function  R.  is  then  found  to  be 


.  ,  ,  /  1  f  L  .  .. 

EEZE-  ,m<8< 

^  m=0  s=0  p=f>  q=C 


(i  +  m)!  Fimp  ®  F<sq  *V  X 


X  (l+pV  ■ 1  H/p(2p_ ()  (p)  [(-1)' +  m'3  uf > _s  cos  +  e;„o>  + 


ipm  iqs' 


+  U“» S  cos  (0.  -  0'  )]  , 

f  '  ipm  fqs'J  ’ 


(25) 


where 


(i  -  2p)  u)  +  mfi  . 


(26) 
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5.  SECULAR  AND  LONG-PERIOD  LINEAR  PERTURBATIONS 


The  Lagrange  planetary  equations  can  be  written  as 


da  _ 

2  8R 

dt  ~ 

na  8M  ’ 

de  _ 

1-e2  8R 

}/l  -e2  8R 

dt 

2  8M 

2  8w  ’ 

na  e 

na  e 

dl  _ 

cot  I 

8R  cosec  I 

dt 

2  *f 

2  du  2,r 

na  wl-e 

na  y  1  - 

dM  _ 

,  2 

1  -e 

8R  _2_  8R 

dt 

n  2 

8e  na  8a  ’ 

80 


na  e 


dw  _ 

cot  I 

dt 

na 

dO 

cosec  I 

dt  2 

na  f  1  - e 

=  MV2 

-3/2  .. 

a  if  it 

M  =  a 

+  f  n  dt  ~ 

3R 


,  f^-e2  8R 


-  31  2  9e  ’ 

2  na  e 


81 


(27) 


-Jte 


8R 

8M 


dt  I  dt 


(28) 


In  the  above  relation,  o  contains  all  perturbations  defined  in  the  fourth  equation 
of  (27).  The  last  term  of  Equation  (28)  contains  only  short-period  terms  and  will  not 
be  considered  in  this  section. 
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Now,  the  function 


l>2 


XT2  2  -I  =■ 

N„  a,  R. 

C  C  < 


(29) 


does  not  depend  on  M  and  is  an  explicit  function  of  time  only  through  r 
considering  I  and  £  constants,  which  is  a  good  approximation. 


c 


v 


C 


,  and  £2^, 


The  integration  of  the  pertinent  equations  can  be  performed  numerically  by  us:,.g 
as  input  lunar  ecliptic  coordinates  -  or,  for  that  matter,  equatorial  coordinates  — 
stored  in  tapes.  This  will  produce  precise  evaluation  of  the  true  lunar  motion.  How¬ 
ever,  such  a  method  can  be  very  expensive  in  time.  A  good  approximation  can  be 
obtained  by  considering  L,,  e  ,  a  ,  and  £  fixed  values  and  M-,  wf,  and  Q  linear 
f>  ictions  of  time,  as  given  in  Section  1,  neglecting  accelerations  of  these  elements. 
AiSO,  an  expansion  in  power  series  of  e^  will  converge  rapidly  owing  to  the  small 
value  of  this  eccentricity. 


Along  these  lines,  we  can  consider  the  expansions 


Gf  .  (e,)  Mfe!  J  , 
?qk'  {'  LcosJV  fqsk/  ’ 


(30) 


where  the  G’s  are  Kaula's  (1961)  coefficients,  which  in  turn  can  be  written  as  Hansen's 
coefficients 


o  =  x-(f+l),  i-2q 

iqk  f -2q  +  k 


(31) 


and 


(32) 
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We  also  remark  that 


0(elk|). 


These  functions  are  given  by 


Gfqk(V"  (1  + 


LJji,'2q+k)et 


j - GO 


f-2q 

f -2q+k,  j 


(33) 


where  the  J.(x)  are  Bessel  functions  with  the  usual  definition 


£  <-»' 


(*/2) 


j  +  2s 


(j  +  S)!  S’ 


s=0 


(34) 


Also,  for  k-j-m  >  0, 


Xf’m  =  /-fl 

k,j  (PC} 


if  +  1  -  m\ 
\k-j-m/ 


F(k-j-f  -1,-m-i  -  1, k-j-mf  l;p  )  (35) 


and,  for  k  -  j  -  m  £  0, 


<-P()'k+i+n’  (-kt/tm)  F(-k  +  i- 


t  -  1,  m-f  -  1,  -k  +  j+  m+  l;Pr) 


(36) 


Here  again,  the  hypergeometric  series  terminate;  that  is,  they  are  polynomials. 
However,  the  G’s  are  infinite  series  that  converge  for  all  e  <  1,  although  for  large 
e.,  the  convergence  is  slow. 

It  follows  that 

(Eq.  eont.  on  next  page) 


lain 

m-0  s-0  p=0  q=0  k=-<» 


£ 


.m 


(-1)  €m«e(f-S)! 


m 


(f  +  m)  ’. 


F  (I)  F  (I  )X 
fmpw  jfsq'  r 
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f  mspqk 


=  (I  -  2p)  u)  +  mh  ±  |{l  -  2q)  o>  +  (i  -  2q+  k)  +  s  , 


dt  =  S*  .  tD 
£  msDak/  j 


'£  mspqk  £  mspqk/  i  mspqk 


In  the  above  relations,  we  are  using  the  secular  rates  of  the  Moon's  motion  as  given 
in  Section  i  and,  for  the  satellite,  as  given  by  the  even  zonal-harmonics  coefficients. 
The  dominant  terms  follow: 


2  2 
I  1-5  cos  I 


w=-4J2\T/  “—2-2 


(1-0 


j2  =  _  2  j  (jef  cos  I 
2  2  W  (l-e2)2 


[,  3  T  /aef  -1  +  3  cos2 1 

M=n[1  +  4J2  Vrj  — ^372-J  * 


The  relation  between  a  and  n  is  the  perturbed  Kepler  lav/ 


22  „  i  3  T 

n  a  =Gm0  1  +  ^  J2 


\2  2  " 
! ]  1-3  cos  I 

'  (l-e2)3/2- 


Now,  let  6.  e.  ,  be  the  linear  long-period  and  secular  perturba- 

’  1  £  mspqk5  ’  1  £  mspqk 

tions  that  are  obtained  +oni  Equations  (27),  integrating  the  right-hand  members 

2  2-1 _ 

(a,  e,I  fixed)  and  substitu  ng  N,  af  jnspqk  for  R'  par^a*  derivatives  enter¬ 
ing  Lagrange's  Equations  (27)  are  given  by 


=  4>.  ,(2p-£) 

jfms.x^k'  r  ' 


[(-!)'  + 


m-s  m, -s  +  m,  s„- 

U£  °£  mspqk  £  £  mspqk  ’ 
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8R 


^geak  =  _*  .  1X1/ + m  - s  uf > -s 

8fi  /  mspqk  £  ’  /  Z 


+  Tjm>  8  s“ 
/mspqk  I  /mspqk 


]■ 

(46) 


3R 


/  mspqk  _  /  ^ 

8a  a  ttZ  mspqk  * 


(47) 


8R 


^S8ES^  =  r  /f  -F1  ) 

81  /mspqk '  /mp  /  mp'  ’ 


and 


(48) 


8R 


/mspqk  _  s-  _  He  -| 

8e  n/mspqk  Ln/p(2p  -/)  /p(2p-/)J  * 


(49) 


where,  if  only  terms  with  positive  powers  are  considered, 


I  =  Y*  (2/.-2i):  2.1 

*/mpW  Lu  i!  (/-i)’.  (/- 


2i  -  2/ 


m-  2i)! 


.  /  -  m  -  2i  - 1  T  v 
sin  1  X 


X  J(Z  -  m  -  2i)  cos2  I  -  j  sin2  i]  cos^  1  I  X 


XE(i'T21%”:i)(-‘)k'q  • 


(50) 


with  the  same  summation  conventions  of  Equation  (9):  for  2p-Z  >  0, 


H 


Zp(2p-Z) 


e  VT7  1 


1  /  +  1 


(2p-£)  p  -  1  +  p2j  H/p(2p-/)  + 


2.) 


+  2(1  +  p2)  (-P)2p  " ' 1  +  1  (2P *  j)  ^  ~  2p2-P|  F(-/ ,  2p  -  2/ ,  2p  -  /  +  l;p2)| 


(51) 
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and,  for  2p  - 1  <  0, 


He 


H 


ip(2p-i) 


-  2(1  +  P2)(-P)i  ‘2p+1(2i  ;^Pp l)  ■Ci/:1jpE2t FH>~2P»i  -2p+2;P2) 


The  following  definitions  are  introduced: 


*,mSpqk  =  <-‘>m  'Sir%r'-  F|sq  V  °iqk<VX 


(52) 


(53) 


*Lspqk=  mspqk  ’ 


(54) 


=  *  /f  — 

i  mspqk  i  mspqk  v  '  ’ 


(55) 


mspqk  ^  ^ 


X +  m-s  T,m, -s  r+  m,  s 

Uf  mspqk  f  i  mspqk  ’ 


(56) 


mspqk 


=  (t  ~2p) 


i  +  m-  s  m,-s  +  ,+  + 

[  '  l  fmspqk^  ( mspqk 


+  mspqk7  f  mspqk  1  ’ 


(57) 


S  mspqk  m 


(-1) 


i  +  m-  s  ^m,  -s  r+ 


,/d; 


i  l  mspqk/  i  mspqk 


,Tm,  s  /j)- 

L(?  ^fmsnak/  f 


'f  £  mspqk'  i  mspqk 


(58) 
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and 


i  mspqk 


=  <-*>' + m ' 8  *7’  '*  + u*m’ 8  iw/i> 


1  mspqk/  l  mspqk 
(59) 


Thus,  we  can  write 


gtmsEgk  (#  =  t 

J  9co 


,0J 


t  mspqk  i  mspqk  * 


JdR 


i  mspqk  r|+_  .  rf2 

9i2  /  mspqk  mspqk  ’ 


/» 


l  mspqk  „  _  f_  g 

9a  a  9 1  mspqk  l  mspqk  ’ 


j3R 


i™m£dt=  <s,c  s 

9e  l  mspqk  i  mspqk  5 


J9R 


i  mspqk  ,  I  s 

91  j?  mspqk  i  mspqk 


The  above  integrals  are  not  valid  if  the  integers 


i-2p,  i-2q+k,  l-2q,  m,  s 


are  simultaneously  zero;  that  is,  we  must  exclude  the  cases 

m  =  0  , 
s  =  0  , 

2p  =  l  -  even  -  2y  , 

2q  =  i  =  even  =  2y  , 

k  =  0  .  (60) 
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They  correspond  to  secular  perturbations  and,  in  this  case, 


Ot>  nfp 

_2v.Q,o,y,y,o  -  ^OAy.y.O  _ 
du  an 

In  order  to  evaluate  the  other  three  integrals,  we  must  consider  (l  =  even=  2y): 


-  21)!  221  ~4y  (2y  -  21) 
(2y  -  i)!  (2y  -  2i)S 


sm 


2v  -2i  -  1  T  oT 

Y  I  cos  I 


(61) 


2y+l 


H, 


2Y,Y,0 


=  F(-2V  -  1,  -2y  -  1,  1;(3  )  = 


(-2Y-1)  (-2Y-1)  2 

- -  (T  ,  (62) 


n! 


n=0 


and 


H 


2Y,  Y>  0 


2  ^  (-2V  -  l)n  (-2Y  -  1) 


efT?  n=l 


(n-  1): 


n  p2n 


(63) 


Therefore, 


3R, 


2YiOA^x£  =  r 

9a  a  2y,  0, 0,  Y,  Y,  0  ’ 


9R, 


2Yl  0j  0,  Yj  Y»  0  _  d  /jr  —  ) 

91  n2Y,  0, 0,  Y,  Y,  0  2Y,  0,  Y  2Y,0,y;  » 


3R, 


gYlg^YxY^O  =  m  /H  _HC  ) 

9e  U2Y,  0, 0,  Y,  Y,  0  n2Y,  Y,  0  n2Y,  Y,  0> 


The  long-period  perturbations  are  given  (excluding  cases  (60))  by 
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. . .  . 


A 


6e  -M2a2-<  T  t  c“ 

1  I  msprk  1  i  \  y  2  9 1  mspqk  *  mspqk 

\  na  c  i 


j  i  ,  =  N2  a2  - - -  (cot  I  C*f  ,  -  cosec  I  )  X 

1  <  mspqk  (  (  2  -/ - 2  '  i  mspqk  i  mspqk  > 

na^  f  1  -  e 


i  mspqk 

> 

x-2  2-i 

1 

9 

'4  \ 

u 

na 

\  6 

N2a2_i 

1 

/  cot  I 

Nc  £ 

0 

w 

na 

Wl-e2 

Pfmspqky  i  mspqk  ’ 


.1  ,  yi+e“  e  g 

?l  mspqk  e  9 1 mspqk/  i  mspqk 


co  _  M2  2-i  cosec  I  „I  ~ 

1  i  mspqk  2  J - 2  ^  mspqk  £  mspqk 

na  fl  -  e 


The  secular  perturbations  are  given  by 


6,  M‘  An  n  =  -n?  a? " 2y  R„  nn  n(H-He)  + 

1  2y,0,0,y,y,°  £  £  na2  *-  6  2Y>°>°>Y>Y>0 


+  4y  R2Yj0,0,Y,Y,0 


3*. 


=  N2  2-2Y  1  L-MJLr; 


j1  w2Y,0,0,y,Y,0  na2  ^  y—2  2y,0,0,y,Y,0 

VT-^e2  -  el 

+  "■■■  -  R0  nn  n(H-He)t  , 

e  2Y,0,0,y,Y,0  I 


(F  -  F1)  + 


(Eq.  cont.  on  next  page) 
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.7?  _  m2  2-2y 

6in2Y,0,0,Y,Y,0  £a<C 


cosec  I  P5 


2y>  0, 0,  y,Y»c 
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6.  SECULAR  AND  LONG-PERIOD  SECOND-ORDER  PERTURBATIONS 


2  2-1 


In  a  second-order  evaluation,  that  is,  if  terms  of  the  order  J0N,a. 

^  C  i 


are  con¬ 
sidered,  it  will  be  necessary  to  take  into  account  the  secular  changes  in  M,  w,  flowing 
to  Jg.  Such  changes  produce  the  largest  higher  order  perturbations  since  they  produce 
amplitudes  that  increase  linearly  with  time.  Let  6^  M,  6^  w,  and  6^  fl  be  these 
perturbations.  Taking  into  account  only  secular  coefficients  and  noting  that  6^  a=  0, 
wt  have 


d  ,c  _  3M  c  —  ,  3M  c  t 
dt  (&2  M)~  ae  6l  e  +  91  6l 1  ’ 


d  .  r  — »  doj  c  ~  3cO  r  J 

dt  (62  ^  “  9e  61  6  91  61 1  ’ 


d  ,  —  9fl  r  —  9fl  r  T 

dt  (52  9e  51  6  91  61  1 


By  considering  Equations  (43),  we  find  that 


d  (62  M)  =  3e  VlTe2  —  ^^  +  3  VT-e2  sin  I  fl  6j  I  , 

1  _  ^  one  T 


dt 


,  r  2t  1 

1-5  cos  I 


^  (62  S)  =  (L  6][  e  +  5  fl  sin  I  6X  I  , 

1  -e 


d_ 

dt 


<s2  n> 


4e 

1-e2 


fl  61  e  -  fl  tan  I  61  I 


» 


where,  again,  w,  fl  are  given  by  Equations  (43).  It  follows  that 


(66) 


Preceding  page  blank 
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6?  M  =  3e  Vl  -e2  1-3  cos  1  u  fs  e  dt  +  3  V*  -e2  n  sin  1  fa  I  dt  , 
1-5  cos^I  ** 

&2  <*>  =  — ^2  w^|*6jedt+5n  sin  I J*  6^  I  dt  , 

62  fi= J&1  5  dt  -  f2  tan!  Js.T  dt  . 


(67) 


Secular  accelerations  do  not  exist,  since  e,  I  have  only  long-period  terms. 
Therefore,  conditions  (60)  have  to  be  excluded.  If  we  consider  the  first  two  Equations 
of  (64),  it  follows  that 


/' 


C"  dt=S" 


'/  mspqk  /  mspqk  * 


/« 


C?  .  dt  =  S?  ,  , 

/mspqk  /mspqk  ’ 


where 


^ mspqk-  L(-1),  +  m_Sur’'6s^spqk/<mspqk'2  + 


+  U/m,8s;mspqk/^mspqk>2 


<*-2p)  , 


3^  = 

/mspqk 


<- 1/ + m  - 8  u“-  -s  s;mSMk/<mspqk>2 + 


+  U™»  8  ST  ,  /(DT  .  )2 
/  /  mspqk/ v  /  mspqk7 


m  . 


(68) 
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*  iiiiiii  iiii 


iMsLnsi>-c 


It  follows  that 


c  —  J  2-1  3y  1  -e2  {  1-ScosT  *[,  2  • 

6o  M,_.=NX  r  7-- - 3-  fl-e  * 


2  fmspqk  C( 


•._£££!_  n  x 


1  -  5  cos  I 


|T7 


oj  2  2-/  3  Q 

X  ^imspqk  Si mspqk  "  \  2  fmspqk  Simspqk  5 

113. 


The  total  long -period  and  secular  perturbations,  including  leading  coupling  terms  with 
J2,  are  finally  obtained  by 

6e  =  6  e  +  6  e  ,  ...  ,  6S7  —  6_£2 

12  1  Z 

Obviously,  the  above  relations  are  not  valid  for  cases  of  critical  inclination  or  satel¬ 
lites  whose  periods  are  commensurable  with  the  rotation  period  (24  h)  of  the  Earth. 
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7.  COMPUTATIONAL  PROCEDURE  FOR  LONG-PERIOD 
AND  SECULAR  PERTURBATIONS 


In  short,  to  obtain  long-period  and  secular  perturbations  due  to  a  term  mspqk 
(Eq.  38),  we  proceed  as  follows: 


-Period  Perturbations 


Compute,  given  means  elements  a,  e,  I,  e^,  1^,  £,  and 

2)  p  (22) 

3>  Fm<v(9) 

4)  Fimp(I)(9) 

5)  Giqk(V  <33)’  (34)’  (35)>  °r  (36) 

6>  HiP(2P-i)^(23)or<24) 

7>  *imspqk<53> 

8)  u®,-,  u-,s  (12a)  Qr  (12b)  frorn  sgcyon  j) 

9)  to,  £2  (43) 

10>  %m  <26> 

11)  M^,  ^  (Section  1) 

12>  %sk  <32> 

13>  CLspqk  <39> 

14>  Sfmspqk  <40> 


Preceding  page  blank 
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15)  n(,  (Section  1) 

16)  w,  M  (43) 

17>  “imspqk  (41> 

W)  FU  (81) 

18)  H/p(2p-()  (51)  OT  (52) 

20>  *'mspqk  (54) 

81)  *|mspqk 
22>  ^mspck*56* 

23)  S(mspqkl59^ 

2“)  Ci“mspqk  <57> 

25>  C?mSpqk  <58) 

‘8)  £„**<“> 

27>  ^/mspqk  <68) 

28)  (element^ ingpqk  (64) 

29)  62  (element^  mgpqk  (69) 

30)  6  (element)  =  6^  (element)  +  (element) 
Complete  long-period  perturbations. 

7. 2  Secular  Perturbations  (i  even  =  2y) 

Given  a,  e,  I,  e(,  I{,  £,  p  (mean  values): 

1)  F2V,0,V  V  (9) 

2>  F2V,0,V(I)<9) 
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3)  G2V;0)0  (e:;>  <33>.  (M)»  <35>,  «  (36) 

4>  h2y,y,o<|3><62) 

5)  U^°(12b) 


6)  t2Y,0,0,v,Y,0(53> 

7)  C2V>  0, 0,  Y,  Y,  0  =  1 

8)  R2y,0,0,y,Y,0(M) 

9>  F2y,0,Y  ®  <61> 

10>  ^Y,V,0  (W  (63) 

“> 

12)  R„..  ofH-H6)  ra> 


2y,0,0,  y,  Y,0 


13)  (element,  secular)  (65) 

Complete  secular  perturbations. 
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8.  SHORT-PERIOD  PERTURBATIONS  OF  LOW  SATELLITES 


During  a  few  revolutions  of  the  satellite,  where  short-period  variations  are  of 
interest,  the  position  of  the  Moon  and  u,  Q  change  little  for  low  satellites.  Then,  as 
Kozai  (1966)  suggested,  we  can  consider  the  elements  of  the  Moon,  w  and  Q,  fixed  when 
performing  the  integrations.  In  this  case,  the  appropriate  expression  for  is  given 
in  Equation  (8),  since  here  it  is  immaterial  what  frame  of  reference  is  being  used  for 
the  coordinates  of  the  Moon.  The  particular  term  R^  will  be  written  as 


t 

- a'  (!)  L,  L,  K cos «'  -2»v  * mt:|  - 

m=0  p=0 


+  S™  sin  [(f  -2p)v  +  mf2]}  , 

where 


(70) 


C 


m 

f 


\Af  > 

f  -  m  even 

»-B  ”, 

i  -  m  rxld  , 

,Rnl 

\Bl  » 

f  -  m  even 

>Am 
l  ’ 

t  -  m  odd  . 

(71) 

The  coefficients  C’s  and  S's  depend  only  on  the  Moon  (5).  The  values  of  r  ,  v\  and 
f2',  given  mean  elements  ar,  e^,  I  .,  £,  and  the  time  T,  can  be  computed  by  consider¬ 
ing  M„,  and  fi  (Section  1),  then  solving  Kepler’s  equation 

'I  1 


sin  Ec  , 


computing  f,  from 


(72) 
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f„  /T+eT  E- 

tan?  =  1/t— ?  tan-? 


computing 


=  a,  (1-e,  cos  E.) 


vc  fc n'  “c  * 


computing  a  ,  6^  from  Equations  (7),  and  finally  a',  6'  from  Equations  (10). 


R  =  ai  f-')  F,  fl)  (C?  cos  0  +  S™  sin  0.  )  , 

imp  Va/  impv;  '  i  ipm  i  ipnr  ’ 


where 


9ipm=  (*-2P)v  +  mn 


The  following  relations  are  easily  established: 


JL2E  =  1  r 
9a  a  imp  ’ 


r\T»  O 

^='^(f)2R,mp  +  ,i(f)R* 


+  -  2P>  *  ©  Fi  mp  »  K"1  sin  %m  + 

+  sr  coseiPm)sinf  (1+f  7^2)  ; 
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8R  t 

“T®  -  a'  (I)  FU«>  (cr  008  %m  +  S/“  si"  %m)  > 


(77) 


where  F*  (I)  is  given  by  Equation  (50); 


8R  i 

-C*’  "  -2«  (I)  «  -  %m  *  ST  -  %J  «  ™ 


9R 


^ma1^  F. 

9S7  Va/  imp 


a)  (-c“ 


sin  0,  +  S®  cos  0.  )  ; 

i  ipm/ 


ipm 


(79) 


and 


9R 


imp_  ,  e  a  „  . 

—  .  -  i  - —  R. _ sin  f  + 


SM 


~2  r  imp 


1-2 


+  ~2P)  a"  ^  (!)  (-Cim  Sin  9fpm  +  *?  tos  0ipm)Fimp(I) 


(80) 


For  short-period  perturbations,  we  make  use  of  Lagrange's  equations  (27)  whert,  R 
is  substituted  by 


2tt 

Rper.  =  B  -  Fi  J  RdM  ■  <8» 

0 

and  integration  is  carried  on  with  respect  to  dt=  dM/n,  considering  all  other  angles 
and  actions  to  be  constant. 

The  computation  of  the  average  of  R  with  respect  to  M  involves  the  following 
integrals: 
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m 


2ir 

I^jP  =  —  J  (f/ cos  (£-2p)fdM=  (l  +  p2)'i_1^Q  29  (P)  , 


(82) 


2it 

r*>P=  —  f  (-)  cos  (i  -2p)f  dM=  (1  +  p2)  ^  xf  q  ’  P  (P)  ,  (83) 

2  2ti  J  w  w 


2ir 

4’P=  2^  J  (F)  C0S  (£  "2p)f  dM=  (1  +  p2)  1  <0M  ^  ((3)  : 

0 


(84) 


2ir 


1^5  P  =  ~  j  ^  sin  ( l  -  2p)  f  sin  f  dM  = 


and 


-  (1  +  P  )  2X0,0  /  ’ 


2tt 


4’P  =  2w  /  (a)  sin  ^  “ 2p* f  sin  f  dM  = 

0 

=  (1  +  P)  \2  X0, 0  2X0,0  )  ’ 


where  X^’  q(P)  is  Hansen's  coefficient  defined  by  Equation  (23)  or  (24). 
It  follows  that,  by  defining 


(85) 


(86) 


and 


Vn=<SinWS<mc0S%m  ’ 

5ipm=  ^  00s%m  +  ST  sin%m  > 


(87) 

(88) 


where  0„  is  given  by  Equation  (26),  we  have 
ipm 
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U  II 

"£/ 


Imp  2tt  J  3M 
0 


dM  =  0  , 


,<2)  .If  S 

imp  2tt  J  3< 
0 


;=<<-2p)a' F  (Dl'j’Ps 


ipm  ’ 


^  n 

■fj 


dM=ma'F<mp<I>!'l’P%m  - 


r(4)  = 


■A/ 


dM-  J  il  J-  F(mpff)  i'.P  0<pm  + 


+  (  S  a<  F(mpW  4’P  Cimp  -  <'  -2p>  a<  F(mpfl)  4’P  C(pn> 


i  1 


(i  2P)a  1  e2  F^mP(I)  IS>  ^P111  ’ 


U  II 

=  1  f 

“2ir  J 


^HdM=2ia^1Fimpa)l{’PCipm  , 


-A  f 

~2rr  J 


^EdM=a'F‘mpa)li’PC( 


Now  it  is  necessary  to  evaluate,  in  closed  form,  integrals  of  the  type 


JR+iJ>  -fgf.*. 

pq  pq  pq  J  Va/ 


for  q  =  0, 1, 2, . . . ,  p  +  1  and  p  =  0, 1,  2, - Introducing  dM  =  (r/a)  dE,  we  e  )tain 
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y 


where 

d)Tfl=  (i-e  e°8  E)^1  , 

eiqf  =  [(cos  E  -  e)  +  177  sin  E]q  , 
7)  ~  |/l  -e2  ,  i  =  \^I  • 


Therefore, 


pq 


-/©’ 


.P+l-Q  _  „,q 


[(cos  E-  e)  +  irj  sin  Ep  dE 


with 


p+ 1 -q  a  0  . 

It  is  easily  found  that 
P+1 

Jpq=  S  y  FWe)(sin  yE  " 1  cos  yE)  +  Kpq°(e)  E  ’  (95) 

Y=-p-l 
(V  *  0) 

where 


=yyv  2-tt(-i)'J-a(p+1"1 

pqy  L-jL-j  L-J  \  k 

a  k  s 


k  W  * 

-  s  /  \Q-k/  \  s 


xe2k+q-a(l^)“'k's  (1--1)8  . 


196) 
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in  which 


a=  l\',  |v|  +  2,  |vi  +  4,  . 
k  =  C ,  1, 2, , 
s=  0,  i,2, ...,  a-k  . 


.,  p  orp+1  , 


Thus, 


JR 

pq 


c  P^l 

a  J  (£)  COS  qf  dM  =  ^ 


Y=“P“1 

(Y*0) 


sin  v  E  +  K  „  E 
y  pqy  Y  pqo 


and 


(97) 


r  pfl 

“/  $ sln  qt  dM  “ '  S 


Y=-p-l 

(Y^O) 


—  K  cos  v  E 
Y  PQY 


(98) 


The  following  integrals  are  then  established: 


A*mp=  ri3^L“EdM=i 

i  -J  ay  2 


a1  — - - F  a) 

imp  v  ' 


vr; 


[(Ji  -l,i  - 2p -  1  “  Ji  -l,|-2p+l)C 


(  TR  -  IR  \  g  1 

\i  -1,1  -2p-  1  i -l,f -2p+l/  ipmj 


+  <'-2 


X  (~Ji  -  2,  i  -  2p  Cf  pm  +  J£  -2,f  -2p  Si*pm)  * 


tR 


(99) 


.  r  aR„ 

^ mp  _ j  imp 


dM=([-2p)a'F 


c._  +j 


R 


^  fmp'  '  V  i,i  -2p  ipm  £,i  ~2p  ipm 


SfDm)  > 


(100) 
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rnamm  ttarrwn m  mur . . 


fpm 


/8R 

-§2£  dM 53  ma'  F/mp(I)  (-J^  .2pC/pm  -  J®,  .2p  V>)  , 


A,mPsj^dM=(a,^l 


7  Imp  ^  Of-  2,  1  -  2p  Clpm  +  J1  -  2, 1  -  2p  Slpm) 


'  i  W*  Of-  1,1  -2p  Clpm  +  ■’J  -  l,l-2p  Slpm)  + 


+  -£aiF,  a) 
2  1  mp  w 


r  ^  ^  c  + 

LV,i-l,i-2p-l  if  -  1,1  - 2p+  1'  ipm 
+  0l  -1,1  -2p-  1  "  Jl-l,l-2p+l^lpm] 


+  5a<7L2(/  -2»  FlmpW  [‘Of,,  -2p-  1  -  JM  -2p+ 1)  Clpm  + 

+  (Jl,i -2p-  1  ~  Ji,i -2p+  1^  Sfpm]  * 

(102) 

Ai5n,P-/^gEdM-^i‘1FflI1pW  (<i-2pWJM-2p%m)  ’ 

(103) 

A6mP  s  J 1 dM =  *  F1  mp  ®  (Jl  l  -  2p  %m  +  J{ ,  i  -  2p  *1  pm)  *  <104> 


Now,  if  we  make  use  of  Lagrange’s  equations  and  consider  Equation  (81),  the 
short-period  perturbations  due  to  a  term  mp  are  given  by 


a  _  XT2  2-1  2  ./mp 

n2aAl 


(Eq,  cont.  on  next  page) 
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7 


[t^?  a'-»  .  (a'"*  .  .f»p  Mj 


Ae 


imp 


n  a  e 


m2  2-i 
N,  a 


AI 


imp 


“7==  cosec  ,  [(Af P  -  if  >p  M)  cos  I  -  (AfP  -  if  >p  M)]  , 


2  2  */!  2 
n  a  f 1 -e 


XT2  2-i 

N- a 


^  mp  '  -  -^Vf-  t1  -  e2>  (ArP  -  Cp  ”)  +  2ae  (A5mP  -  fmp  M)1  +  A'M 


mp 


imp 


i  mp  f 


m2  2-i 
N„  a 


Au 


imp 


[-  e  cot  I  (A'”*  -  if  »p  m)+  ,1  -e2,  (a'”»  -  if  >p  m)]  , 


2  2  *[.  2 
n  a  e  F 1  -e 


imp 


Aft 


imp 


nJ  a^_i  cosec  I 

i  i _ 


2  2  i/i  2 

n  a  fl-e 


( AimP_  T  (6)  M) 
VA6  Aimp  '  ’ 


(105) 


where 


A'M 


imp 


om2  2-i 

3Nx!l  Ja< 


mp 


2  2 
n  a 


dM 


(106) 


remains  to  be  evaluated.  This  last  involves  the  evaluation  of  the  integral 


L  =  |j  dM=  jj  (1  -e  cos  E)  dE  =  +  i  L*  . 

pq  J  pq  J  pq v  '  pq  pq 


It  is  readily  found  that 


P+1 


L“  =  (V1  dM=  -  Y 

pq  J  pq  L-J 


Y=-p-l 

(Y^O) 


7  KpqV  C°S  V  E  +  V  T 


(Eq.  cont.  on  next  page) 
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r  - 

e  Km0  (E  sto  £  +  cos  B)  + 1  ]T  ^71)  cos  (Y  + 1)  E  + 


Y=-p-l 

(y*0,y*-i) 


r '  - 

I  E  wn)Kpqvcos(Y-,)E  • 


Y=-p-l 

(Y^OjY*1) 


Lpq=  JJwdM  = 


E  7  Kpqv  Sia  vE  +f  ^pql  •  "p.q,-!'®  + 

Y=-p-l  Y 
(Y^O) 


pfl  P+1 

1  V*  ■■;  1  K  sin(Y  +  l)E+|  Y  -~-rr  K  sin(Y-l)E 

2  Z-rf  y(y  +  1)  WY  2  y(Y"1)  WY 

Y=-p-l  Y=“P“1 

(Y*0,y*-1)  (y^0,y^!) 


Therefore, 


A/„  3  „2  2-1  la1'2  e  „  m  [7T I 

AMimp  "2NCaC  Q2  y— ^mp^  L^Li -l,i-2p-l 


-  l,f  -2p+  1 


)Cipm  +  (Lf  -l,f  -2p-l  "  Li-l,l  -2p+l)Sfpm]  “ 


(I-2P)  4^-  VT7  F/mp(I)  (l? . 2> ,  .  jp  S, 


i  -2,f  -2p  ipm 


C.  )  ■ 
fnrrv 


c--1*-  v --  i— ■ .  v- 


9.  COMPUTATIONAL  PROCEDURE  FOR  SHORT-PERIOD  PERTURBATIONS 

OF  LOW  SATELLITES 

Corresponding  to  a  particular  term  R^mp  (Eq.  74),  the  short -period  perturbations 
(l,  oj,  Q  fixed)  are  computed  as  follows. 

Given  mean  elements  a,  e,  I,  e^,  1^,  and  E,  we  compute 

iv  xt2  2-i 

1)  N.  a. 

'  C  C 

2) p  (22) 

3)  Given  compute  E^,  f 

4)  Given  co^,  and  a^,  compute  r^,  v^ 

5)  a(>  6C  (7) 

6)  a',  6'  (10) 

7)  P™  (sin  6')  c?s  m  6 ' 

8)  A“  B™  (5) 

9)  C™,  s“  (71) 

10)  Given  M,  w,  and  Q,  compute  E 
%m<26> 

12>  Cipm,  *ipm  (87),  (88) 

13>  Fimp^  <9> 

!4)  Fjmp(I)  (50) 

15)  X^’j0((3)  (23)  or  (24) 


k 


j 


£ 

£  -2 
I  - 1 
£  -1 
1-1 
£ 
i 


I  -2p 
£  -2p 
£  -2p 
I  -2p -1 
i-2p+i 
i  -2p - 1 
l  -2p+  1 


16)  KmY (e)  (96)>  v  =  -p-i,  -p,  •  •  * ,  p,  p+  1 


p  q 


l  - 1 

£  -2p - 1 

£  -  1 

£  -2p+l 

£  -2 

£  -  2p 

£ 

£  -  2p 

£-\ 

£  -2p 

£ 

£  -2p-  1 

£ 

£  -2p+  1 

17)  iJ’P  4,P,  4’P  (82)  through  (86) 

18>  fmp>  Cp>  Cp  (83)  «toough  (93) 

19)  J11  ,  J1  (97),  (98)  (same  range  for  p,q  as  in  (15)) 

pq  pq 

20)  A^mp,  A*mp,  . . . ,  Apmp  (99)  through  (104) 

21)  Lpq’  Lpq  (107)’  (108) 

p  q 


£  -  1 

£  -2p-  1 

£  -1 

£  -  2p  -f  1 

£  -2 

£  -2p 

44 


22> 

23>  >  A‘W105> 

Short -period  perturbations  completed. 
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10.  SHORT -PERIOD  PERTURBATIONS  OF  HIGH  SATELLITES 
WITH  SMALL  ECCENTRICITY 


When  the  satellite  is  high,  for  example,  close  to  a  24-h  period,  the  Moon  can  no 
longer  be  considered  fixed  during  a  few  revolutions  of  the  satellite.  Here  we  consider 
also  the  variations  of  to,  O,  in  contrast  to  what  we  have  done  in  Section  8.  In  this 
case,  the  integrals  found  in  that  section  have  to  take  these  variations  into  account. 
This  can  be  done  only  if  the  eccentricity  of  the  satellite  is  small  so  that  power  series 
in  e  will  converge  rapidly.  Thus,  the  disturbing  function  obtained  in  Equation  (19) 
has  to  be  developed  in  terms  of  M,  M^ . 


The  following  expansions  are  well  known: 


E  Glqk<eC>  cos  0Lqk 

k=-OG 


(110) 


and 


sin 

cos 


[{i  -  2p)  v  +  mfi] 


>WC) 


sin 

cos 


impj 


(111) 


where 


£  sqk 


(i  -2q)  ^  +  (f  - 2q  +  k)  +  s(^  +  ~) 


(112) 


eimp.=  (i -2p)w  +  (^ -2P  +  j)M  +  , 


(113) 


and  H.  .  (e)  are  Kaula's  (1962)  coefficients.  These  can  also  be  written  in  terms  ol 

r  J 

Hansen's  coefficients  by 


Preceding  page  blank 
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The  coefficients  (e^)  have  been  defined  in  Equations  (33),  (34),  and  (35)  or  (36). 

The  classical  expressions  for  Kaula's  coefficients  are  (e.g,  Plummer,  1960,  p.  45): 


Hipj(e)=(l  +  p2)'i"1  ^  Ji[(i-2p+3)e]xJ>_i2"p2+Pjji  =  0(elJi)  ,  (114) 

i=-oo 

where  J^(x)  are  the  usual  Bessel  functions  and  xj^  ^(P)  are  given  in  terms  of  hyper¬ 
geometric  series  (which  always  terminate),  as  follows:  for  k-i-m  >  u, 


<-p> 


m 


-  m  +  1;P2) 


(115) 


and,  for  k  -  i  -  m  :£  0, 


(-P) 


-k+i  +  m  (  l  +  l  +  m 


(-k+i  +  m)F(~k+i“i  _1»  -k+i+m+l;P  ) 


(116) 


It  follows  that 

till  00 

IIL  E 

m=0  s=0  p=0  q=0  k=-oo  j=-oo 


R 


fmspqkj 


where 

R<mspqkj~  ((  +m)!  a  FJmpl1*  F/sq  X 


(Eq.  cont.  on  next  page) 
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X  Gf  qk  <e«>  H<p)  ‘e>  [<-'»'  +  3  «?’  'S  cos  mpj  +  6/  «,k> 

+  ufm’scoS(e<mpre«aqk)]  . 


In  what  follows,  it  will  be  clear  that  we  should  have 


l  -  2n  +  j  *  0 


We  can  make  use  of  the  fact  that 


=  — F(a+1,  b+1,  c  +  1 ; x) 


*i<x>  2 


=  2  [Ji-i<x>-Ji+i<x»  » 


tnd  arrive  at  the  following  relations: 


,yYl,  m 


i<  k-m 


e  ?1  -e  k 


l[Wke>--WHXM 


+  E  |ji<te>zM,n+2p-l[Ji-i<lie>-Ji+i<keilxi;r!  - 


isk-m 


where,  for  i<  x  -  m, 


49 


ax?»m 


k-i-m+  1  (l  +  l-m\ 
\k-i--m/ 


X  (k“i“k 1}  r<k-<  k-i-m+2;p2) 


(119) 


and,  for  i  >  k  -  m, 


ax?’m 

<:?  *  -  -<-k+‘+-»K:r  -  ^)-k+i+m+i  c^rjx 


x 


(-k+i-i  -  l)(m-i-  1) 
-k+  i+  m  +  1 


F(-k+i-f,  m-i,  -k  +  i  +  m  +  2;P  ) 


(120) 


Finally, 


3H(  , 
3e 


p 

e  |/l  -  e2 


vM-2\) 

i-Zp+j  • 


(121) 


Let  us  consider  the  definitions 


Djf mspqkj  =  ^  -2p)«  +  (*  -2p  + j)  M  +  mf2  +  (i  -2q)  +  (f  ~2q  +  k)  +  sf^ 

(122) 

and 

D< mspqkj  =  "2p)  “  +  (i  -2p+^  M  +  mr2  -  (f  -2q)  -  (f  -2q+k)  . 

(123) 

We  easily  establish  that 


«(!)  -  f9Ri  mspqkj 

f  mspqkj  J  3M 


=  R 


W 


ir.,  s 


1  mspqkj  i^m, 


u“’S(l  -2p  +  i)/D‘mspqkj 

-  u“’  's  (i  -  2p  +  i)/DJmapq|tjS  ’ 


(124) 
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B 


(2) 

£  mspqkj 


s  fBR|  mspqkj  dt  = 

\j  J  3uj 


\u(m’  S  -  u“>  s  ,/  -2p)/Djmspqkj  / 

mswkj|um>-s_ur)-S(,.2p)/D,ms^ 


(125 


B 


(3) 

£  mspqkj 


=  f  ^mspqkj  = 
ij  ~)  BCl 


=  R 


iT,m,  s  TTm,  s  * 

’  £  £  7  £  mspqkj  f 

£  mspqkj  i  m,-s  m,-s  /  +  i 

U£  U£  /L>1  mspqkj  i 


(126) 


R(4)  .  f9R£  mspqkj  . ,  a  (5)  r„  _  N  l 

B£ mspqkj  J  de  £  B£mspqkj^  £pj  Ipj^  ’ 


(127) 


and 


b(5)  3  f  dt  = 

B£  mspqkj  J  3a 


TTm, -s  m,  -s/n+  < 

i  \Lf  L£  'D£  mspqkj  j 

a  Rf  mspqkj  ^  u“»  8  -  U™>  S/D~ mspqkj  ^ 


cos  —  sin 


B 


(6) 

£  mspqkj 


3R, 


(ulvl  mspqkj  ,u_  a  p(5)  iF  _Fl  \ 
i  “J  31  dt  £  B£  mspqkj  'B£mp  Imp* 


(128) 


(129) 


A'M.  ,.=  --7  N,  a. 
£  mspqkj  2 


3  .t2„2-£  JB(1)  i.dt  = 


£  mspqkj 


3  m2  2-f 
~2  N,  a  ,  R 


m,  s  _  m,  s  / 


H  i  TT^>  -S  ~S  //n+  i 

U*  '  '  £  mspqkj  \ 

cos  —  sin 


£  mspqkj  iU£  ’  —  U£ 

*  s\r\c>  __  oi  n 


(130) 
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Finally,  the  short-period  perturbations  are  given,  for  all  terms  for  which 


i  -  2p  +  j  #  0  , 


(131) 


by 


2  2-1  2_  R(l) 

j?  mspqkj  \  na  imspqkj  ’ 


Ae 


imspqkj  C  i 


AI _ ,,  =  Nt  at 


am  , .  =  -n;  a; 
i  mspqkj  C  1 


AuJi  mspqkj  Nif  \ 


i  ^1-e2 

V; 

2 

t  1 

na  e 

1  cosec  I 

i 

2*/,  2 
na  f  1  - e 

j 

-i  1 

2 

(i- 

na  e 

i  1 

2  •/T” 

2 

na  e  y  1  - 

■e 

i  cosec  I 

2  n(l)  _  n(2) 


—  B' 

imspqkj  impsqkjl  ’ 


j(2)  cog  j  _  g(3) 

3 i  mspqkj  C  Ui  mspqkj 1  ’ 


i  mspqkj  -1  i  mspqkj 


+  A'M 


i  mspqkj  ’ 


-e  B(.6^  ,  .  cot  I  +  (1  -e2)  B^4)  .  . ,  , 

imspqkj  '  ’  imspqkj1  ’ 


.  „  _  «  ■*-  vvww  a.  y-j  / 

Ai“i mspqkj  “  \  2  / - g  imspqkj  ’ 

na  V 1  - e 


(132) 


which  completes  the  calculations. 
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11.  COMPUTATIONAL  PROCEDURE  FOR  SHORT-PERIOD  PERTURBATIONS 
OF  HIGH  SATELLITES  WITH  SMALL  ECCENTRICITY 


The  sequence  of  calculations  to  obtain  short -period  perturbations  due  to  a  particular 
term,  R^mspq}cjj  °f  the  disturbing  function  (see  (117))  is  now  given. 


Given  mean  elements  a,  e,  I,  M,,  M,  cj,  O,  £,  Mf,  wf,  £2,, 

P,  compute,  for  any  term  (i ,  m,  s,  p,  q,  k,  j),  1  -  2p  +  j  *  0: 


M. 


j,  U, 


"  4qk’  e£mpj  <112>’  <113> 

2)  J.  [(i  -2p+ j)e]  (Bessel  function  (34))  to  the  approximation  required 

3)  (115)  or  (116) 

4)  (9) 

5>  W«)(9) 

6)  Gfqk(eC)  <33>’  <34>’  (35)j  aild  (36) 

7)  Hipj(e)(114) 

8)  Y^’  (P)  (119) 

9)  z£’™(P)  (120) 

10)  m(P)  (118) 

11)  Nfpj(P)(121) 

12>  Df  mspqkj  <122>’  <123> 
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13)  U™» S,  ujn»  "S  (12) 

14)  Bimspqkj»  i=  1»2»-“’6  (125)  through  (129) 

15>  A'Mimspqkj  <130> 

16)  mspqkj’  *  *  * »  ^imspqkj  (132) 

Complete  short -period  perturbations. 
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12.  REMARKS  ON  SOLAR  PERTURBATIONS 


The  previous  formulations  apply  as  well  to  solar  perturbations,  which  are  about 
of  the  same  order  of  magnitude.  In  fact,  for  the  Sun, 


R  = 


(COS  $q)  , 


(133) 


so  that 

mO  _2  „3  J2  „3 
GmO  =  m.  4-  mQ  +ln0  °0  aO  nO  O 

or 

Gm0=N^a^  ,  N2  -  0.  75  X  10~5  rev2  day'2  , 

which  is  of  the  same  size  as  N2.  For  the  Sun,  we  have  (to  the  mean  equinox  of  date): 

u  =  281°13’15V0  +  6 189V 03  T  +  1V63  T2  +  0V012  T3  , 

MQ=  358°28'33V0  +  129596579V 10  T  -  0V54  T2  -  0V012  T3  , 
eQ  --  0.  01675104  (supposed  constant)  , 
aQ  =  1.00000129  (astronomical  units)  , 

Hq  =  3548V 19283  day"1  . 
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We  can  consider  Iq,  £2q  to  be  zero.  The  mean  inclination  with  respect  to  the 
equator  is  £.  For  that  matter,  it  cou’d  be  considered  a  function  of  time,  but  such 
precision  is  hardly  necessary.  The  disturbing  function  is  given  by  Equation  (8),  while 
(2)  is  written 

h=£44-S  ■  <134> 

is  2 

The  transformation  (10)  is  not  necessary,  so  that  the  coefficients  A™,  B™  (Eq.  5) 
retain  their  original  form  by  using  I'  =  £,  the  inclination  of  the  orbit  of  the  Sun  with 
respect  to  the  equator.  It  follows  that 

+ 1  i  i  i 

£  'n,  ^F/mp<I>IW£>X 

m=0  p=0  q=0 

X  cos  [(i  -2p)v  -  (i  -2q)vQ+  mil]  .  (135) 

The  secular  and  long-period  part  of  this  is 

Va<  E  E  L  ‘m  <7^Fin,p<I>,W£)X 

m=0  p=0  q=0  k=-oo 

X  cos  [(i  -  2p)  w  -  (i  -  2q)  WQ  -  (i  -  2q  +  k)  +  mil]  ,  (136) 

which  for  the  Sun  is  used  in  place  of  Equation  (37).  More  precisely: 
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R/mpqk=  ‘m  §7i$  V  W®  '  ‘  Hip(2p-«)  »>X 

X  ^qk^O*  Cipqmk  »  (137) 

where 

Cfpqmk  =  COS  ^  _2P)  u  “  (*  -2q)“0-  (*  -2q  +  k)MQ  +  mfl]  .  (138) 

From  this  point  on,  all  formulas  developed  for  the  Moon  can  be  easily  adapted,  a  task 
not  worth  undertaking  here.  The  complete  expressions  are  given  by  Kaula  (1962),  and 
the  computational  procedure  is  similar  to  the  ones  given  in  the  previous  sections. 
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